Lattice Determination of Heavy-Light Decay Constants 
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We report on the MILC collaboration's calculation of /a, fs,, fo, fo^ and their ratios. Our central 
values come from the quenched approximation, but the quenching error is estimated from Nf ~ 2 
dynamical staggered lattices. We use Wilson light valence quarks and Wilson and static heavy 
quarks. We find, for example, /s = 157 ± 11 tf MeV, /b7/b = 1.11 ± 0.02 to.oa ± 0.03, 
fo, = 210 ± 9 tf MeV and /s//o, = 0.75 ± 0.03 t^ l^ iaoo. where the errors are statistical, 
systematic (within the quenched approximation), and systematic (of quenching), respectively. 



PACS numbers; 12.15.Hh, 12.38.Gc, 13.20.-v 

The neutral B meson (Bd), a bound state of a d quark 
and an anti-6 quark, is known to mix with its antiparti- 
cle, Bd- In the Standard Model, Xd, the ratio of the mix- 
ing oscillation to the decay rate, is proportional to the 
absolute square of the fundamental quantity Vtd- How- 
ever, despite the fact that the Xd is well measured [|], 
Vtd remains poorly determined because the proportion- 
ality constant between Xd and \Vtd\'^ depends on non- 
perturbative strong interaction effects. These effects are 
parameterized by fs , the pseudoscalar decay constant of 
the Bd meson, and Bb, the corresponding "bag param- 
eter." Accurate computations of Jb and Bb therefore 
put tight constraints on the Standard Model. Similarly, 
a measurement of for Bg mesons would determine a 
second fundamental quantity, Vts, if fs^ and Bb^ were 
known, or \Vtd/Vts\, if the ratios /bJ/b and BbJBb 
were known. 

Lattice QCD offers a way to compute quantities like 
/b and Bb from first principles. Here, we present a com- 
putation by the MILC collaboration of the decay con- 
stants /s, fB,, fo, fos, and their ratios. Ref. j|] gives 
additional details; preliminary results were described in 
Refs. 

TableT^ shows the lattice parameters used. See 
for details of the lattice generation, gauge fix- 
ing, and determination of the quark propagators. We 
compute "smeared-local" and "smeared-smeared" pseu- 
doscalar meson propagators in each of three cases: heavy- 
light, static- light, and light-light (with degenerate masses 
only). Light Wilson quark propagators are computed 
by a minimal residual algorithm for three values of the 
hopping parameter, giving light quark masses (wg) in 
the range 0.7ms<mq<2.0ms, where nis is the strange 



quark mass. The light-light pseudoscalars are used to set 
the scale (through /^r) and to find the physical values of 
Kud and Ks, the hopping parameters of the up/down and 
strange quarks. We determine Ks by adjusting the de- 
generate pseudoscalar mass to ^J2m?^ — rn^, the lowest 
order chiral perturbation theory value. We also com- 
pute smeared-local light-light vector meson propagators, 
which we use for alternative determinations of the scale 
(through rUp) and Kg (through m^). 

Heavy quark propagators are computed by the hop- 
ping parameter expansion |^ . Because of practical limi- 
tations to this approach |^ , we sum the sink point of the 
smeared-smeared correlators only over a subset of points 
in a spatial volume. This means that intermediate states 
of nonzero 3-momenta can contribute. For the heavy- 
light mesons studied here, these higher momentum states 
are suppressed sufficiently at asymptotic Euclidean time 
t by their higher energy, although in the largest volumes 
(sets N and O) this can require tmin/a as large as 25 (a 
is the lattice spacing). 

The static-light mesons have no such suppression. 
However on our smallest volumes (sets A, C, D, F, G, H) 
the contamination by higher momentum states is small 
(ss 0.7%, which we estimate using static-light wavefunc- 
tions from Ref. ||^). On all other sets the contamination 
is expected to be large. We therefore have performed a 
dedicated static-light computation on those lattices, with 
relative smearing functions taken from |Q and zero mo- 
mentum intermediate states enforced by a complete EFT 
sum over spatial slices. In addition, the dedicated static 
light computation has been run on set A (because the 
plateaus from the hopping method proved to be poor) 
and set G (as a check of the hopping method). On the 
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latter set, the two methods give consistent results. 

For all pseudoscalars, we fit the smeared-local and 
smeared-smeared correlators simultaneously and covari- 
antly to single exponential forms, with the same mass in 
both channels. We vary the fit range (in t) in each chan- 
nel over several choices that have reasonable confidence 
level (CL). Combining such choices for the light-light, 
heavy-light and static-light cases, we have approximately 
25 different versions of the analysis on each data set. Our 
central values are taken from the version which has the 
best blend of high CL and small statistical errors. We 
then find the standard deviation of the result over the 
other versions and add it in quadrature with the raw 
jackknife error of the central value. The resulting error 
will be called, henceforth, "the statistical error." 

We employ the EKM norm |^ throughout. In the 
heavy-light case we also adjust the measured meson 
pole mass upward by the difference between the heavy 
quark kinematic mass (7712) and the heavy quark pole 
mass (mi) as calculated in the tadpole- improved tree 
approximation fixing the mean link from Kc- We 
use the one-loop tadpole-improved, mass-dependent per- 
turbative renormalization of the axial current with 
coupling Q!v(3.4018/a) defined in terms of the plaquette 
pof . We adjust the result of |^ for our matching point 
(to2 rather than mi) and for our choice of the mean link. 
Our central values use "scale choice i": q}f^ = 2.32/ a for 
the heavy-light corrections and qj^ — 2.18/a for the 
static-light corrections The heavy- light scale was 

calculated in the masslcss limit; however, since it dif- 
fers little from the static-light scale, it seems reasonable 
to use it for all mass values. The effects of two other 
choices of scale (it q*fj]^ — q^j^ — 1/a; in: q*^^ — 4.63/a, 
q*g^ — 4.36/a) give an estimate of the perturbative errors. 

In our chiral fits we use am2 as the independent vari- 
able rather than the more standard Although the 
two are formally equivalent at this order in a, m2 has the 
advantage that CL of linear fits to Mqq and fqq are quite 
good. (Here, Q is a generic — possibly static — heavy 
quark, and q is a generic light quark.) Further, linear fits 
to also have reasonable CL for quenched (3 > 6.0. For 
771^, however, the CL of the linear fits is uniformly very 
poor whether 1/k or am2 is the independent variable. To 
study this problem in more detail, we have examined the 
pseudoscalar mesons for 6 light quark masses at /3 = 5.7 
on additional lattices (set "5.7-large") of size 12^ x 48 
(403 configurations), 16^ x 48 (390 confs.), 20^ x 48 (200 
confs.), and 24^ x 48 (184 confs.). The lightest two me- 
son masses in this set (« 385 and w 515 MeV), are below 
those used in the full computation. On set 5.7-large, lin- 
ear fits of m^ vs. either 1/k or am2 are still poor, but 
quadratic fits are good. Indeed, a quadratic fit of m^ 
vs. am2 using the 5 heaviest masses goes right through 
the lightest mj in all but the 12'^ volume. For our cen- 
tral values we thus employ quadratic fits vs. am2 for mj , 
and linear fits vs. am2 for Mg^, and fqq. We call 



this "chiral fit /." Three other fit choices (//: all linear; 
///: mj and /^r quadratic, all others linear; IV: mj, /tt 
and fqq quadratic, all others linear) are used to assess 
the systematic error. One set, F, undergoes very large 
(^ 50%) variation when the chiral fit choice is changed, 
possibly because of finite size effects. Set F is therefore 
dropped from further analysis. 

To find Jb on a given data set, we divide out the per- 
turbative logarithms Q from fqq ^Mqq, fit to a polyno- 
mial in 1/Mqq, interpolate to ms, and then replace the 
logarithms. We do three versions of the polynomial fit: 
(1) a quadratic fit to the mesons in the approximate mass 
range 2 to 4 GeV ("heavier heavies") (2) a quadratic fit 
to the mesons in the approximate mass range 1.25 to 2 
GeV ("fighter heavies") (3) a cubic fit to the mesons in 
the approximate mass range 1.25 to 4 GeV. We include 
the static-light point in all three fits. We use range (1) 
in central values for fs and /sj range (2), for and 
fu^ . The alternative ranges go into the systematic error 
estimates. 

The final extrapolation is in lattice spacing. Since the 
Wilson action's leading errors are 0(a), we attempt a lin- 
ear extrapolation in a for all our quenched results. Figure 
|l| shows the extrapolation for /s, with the central choices 
of the perturbative scale (choice i) and of the chiral fits 
(fit I) . An alternative possibility, with which the data are 
also consistent, is that the 0{a) effects are small enough 
for 6/52 > 6.0 (a<0.5 GeV"^) that one may extrapolate 
with a constant fit in this region. For the decay con- 
stants, both fits have acceptable CL, but the constant fit 
is better. However, for and fos/fo, the linear 

fits (CL « 0.6) are much better than the constant fits 
(CL « 0.1). (See Ref. for plots of the ratios and ad- 
ditional details.) Since it would be inconsistent to treat 
the decay constants as independent of a, yet fit the ratios 
linearly, and since we in any case expect significant 0{a) 
errors for Wilson fermions, we take the linear fits to the 
quenched results for our central values. The differences 
with the constant fits are included in the systematic er- 
rors. At this point, the dynamical Np — 2 data is not 
good enough to extrapolate to the continuum, even for 
unphysically large dynamical quark mass. We use the dy- 
namical data only to assess the error due to quenching. 
The systematic errors are computed as follows: 
(1) The three largest sources of error within the 
quenched approximation are the continuum extrapola- 
tion, the chiral extrapolation, and the 0{a'^) perturba- 
tive corrections (as estimated from a change in scale in 
the 0{as) terms). With our data, these errors cannot 
be computed independently. For example, when the chi- 
ral extrapolations are changed to fit IV (see Fig. ^) , the 
difference between the linear and constant (not shown) 
continuum extrapolations gets smaller (15 MeV instead 
of 23 MeV). Further, while the systematic error in the 
final results would be very small if the only source of 
uncertainty were the next higher order perturbative cor- 
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rection (©(a^)), this is not the case once the interplay 
between perturbative uncertainties and other continuum 
extrapolation errors is included. Indeed, assume there 
exists "perfect" data which are linear in a with slope and 
intercept as in Fig. |l|, and then add on an ©(a^) correc- 
tion with a coefficient chosen to give the same change at 
(3 = 5.7 as would be produced by reducing and qj^ 
to 1/a (choice ii). Although this gives a 17% change at 
(3 = 5.7, a linear extrapolation of the changes to a = 
results in a residual error of less than 1%. With the real 
data, however, changing q* to choice ii raises the linearly 
extrapolated value by 10% (see Fig. ||); while it reduces 
the constant fit by 3%. Choice Hi reduces the linear value 
and raises the constant value by ~2%. We therefore es- 
timate the errors from the continuum extrapolation, chi- 
ral extrapolation and perturbative corrections together. 
We compute each quantity 24 times (2 continuum ex- 
trapolations X 4 chiral fits x 3 scale choices), giving a 
central value and 23 alternatives. The alternatives are 
divided into two groups depending on whether the result 
is greater or less than the central value, and the standard 
deviation of each group about the central value is then 
taken as the positive or negative combined error. 

(2) The "magnetic mass" rria, which divides the chro- 
momagnetic interaction in the effective non-relativistic 
Hamiltonian for Wilson fermions, is not equal to the ki- 
netic mass 7712 @- This introduces an error at fixed a of 

0{{Cmag " 1) AgcD /Mg,) , where Cmag = TO2/m3. The 

error is not completely removed by the linear extrapola- 
tion to a = 0. Following [Q, we estimate the residual 
error by using the tree level expression for Cmag (with 
our values of 07712) and extrapolating Cmag linearly in 
a. With our preferred choices for the mass range in the 
Jqq ^Mqq fit, this gives an error of ~ 2% for /s and 
^ 3% for /d. The error on /s can be reduced to less 
than 1% by switching to the "lighter heavies" (-)- static) 
mass range: the static-light point, for which 7773 7^ m2 
is not an issue, becomes particularly important in this 
case. In practice, we assess the errors due to 7773 ^ 7772 
as the larger of: (a) the 2 or 3% model estimate with our 
preferred mass ranges and (b) the actual difference in the 
final result caused by switching from "heavier heavies" to 
"lighter heavies" or vice versa. 

(3) Our preferred fits of fgq ^Mq^ vs. l/Mqq are 
truncated at quadratic order. A scale of ~0.75 GeV for 
l/Mqq is expected in the omitted cubic term, since this 
is roughly the scale size found in the linear and quadratic 
terms. We calculate that the existence of such a cubic 
term in the data would lead to an error, in the analysis 
that uses only quadratic fits, of ~ 1% in the decay con- 
stants. In practice we estimate this error by changing to 
cubic fits (using the entire mass range 1.25 to 4 GeV); 
the errors found are indeed <1%. 

(4) The finite volume effects are estimated by com- 
paring results on sets A (spatial size ^1.2 fm) and B 
(~ 2.5 fm) and applying the fractional difference to the 



final results. Set A is smaller than all other quenched 
lattices; B, much larger. Therefore the difference should 
give a conservative bound on the finite volume error. In 
practice, we take the larger of: (a) the difference when 
all quantities are computed individually on sets A and B 
and (b) the difference when all light-light quantities are 
taken from set 5.7-large. Since there is some cancellation 
of error between fqq and /tt, (b) is generally larger. We 
find an error of ~ 2-3% on decay constants, ~ 4% on 
and ~l-2% on other ratios. 

Errors (2)-(4) do not have definite signs and appear to 
be largely independent of each other and of error (1). We 
thus take the error within the quenched approximation 
to be the sum, in quadrature, of errors (1) through (4). 
For decay constants, error (1) always dominates; while 
for the ratios, error (2) (and for fs/fDsi (4)) is (are) 
comparable to (1). 

(5) The quenching error is estimated in three ways: 
(a) We set the scale by using 77ip instead of (b) For 
quantities involving the strange quark, we fix from 777^ 
instead of the pseudoscalars. (c) We compare the results 
from the weakest coupling Np = 2 lattices (sets G and R) 
with the quenched results interpolated {via the linear fit) 
to the same value of the lattice spacing (a w 0.45 GeV~^ 
— see Fig. |]). We average this difference over 12 analysis 
choices (4 chiral fits x 3 scale choices), plus (for strange 
quark quantities) the preferred choices but with Kg fixed 
from 7770. For the decay constants, this difference has a 
definite sign over all 12 or 13 choices. We then take the 
signed error (c) to be just the average difference. How- 
ever, for some of the ratios, the standard deviation of 
the difference is larger than the average difference. In 
that case, the positive (negative) error (c) is taken to be 
average difference plus (minus) the standard deviation. 
Finally, the quenching error in the positive or negative 
direction is defined to be the largest of errors (a), (b), and 
(c) in that direction. In almost all cases, (c) is largest. 

Note that our quenching error estimate is still rather 
crude. Our Np = 2 simulations are not "full QCD" be- 
cause they are not extrapolated to the continuum or to 
the physical quark mass, and they do not have a dynam- 
ical strange quark. For these reasons we prefer to quote 
the central values as the quenched results and to treat 
the difference (c) as a signed error, not a correction. (See 
1^ for further discussion.) Additional dynamical simula- 
tions, which we hope will allow us to improve this situ- 
ation, are in progress. Note, however, that the sign we 
find of the difference is what is expected from intuitive 
arguments about the wave function at the origin ]l^ . 

We then have: 

fs = 157 ± 11 tf MeV 
/b, = 171 ± 10 If tf MeV 
/z5 = 192 ±11 if IJ^ MeV 
fo, =210 ±9 tf MeV 
/bJ/s^I.II ± 0.02 i°:g^ ± 0.03 
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foJfD = 1.10 

fs/fD^ = 0.75 
IbJ/d^ =0.85 ±0.03 



I n 09 +0-04: +0.02 
=1= U.UZ _o.o2 -0.03 

-I- n n"? +0-04 +0.08 

in u.uo _o.02 -0.00 
+0.05 +0.05 
-0.03 -0.00 ■ 



(1) 



where the errors are statistical, systematic (within the 
quenched approximation), and systematic (of quench- 
ing), respectively. The result for fo, is consistent with 
the experimental value ^ of 241 ± 21 ± 30 MeV. Our 
quenched approximation values are consistent with re- 
cent quenched results using improved actions [p^ , p^ . 
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FIG. 1. /b vs. a for preferred choices: chiral fit /, pertur- 
bative scale i, and "heavier heavies" mass range. The scale 
is set by The linear fit to all quenched points (solid line) 
gives the central value. 
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FIG. 2. Same as Fig. ^j, but quenched results only, with 
alternative analysis choices. Some points have been displaced 
slightly horizontally for clarity. 
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TABLE I. Lattice parameters. Sets F, G, and L-R 
use variable-mass Wilson valence quarks and two flavors of 
fixed-mass staggered dynamical fermions; all other runs use 
quenched Wilson quarks. 



set 


(3 (aniq) 


size 


# confs. 


A 


5.7 


8^ X 48 


200 


B 


5.7 


16^ X 48 


100 


E 




12-'' X 18 


100 


C 


6.0 


16^ X 48 


100 


D 


6.3 


24^ X 80 


100 


H 


6.52 


32^ X 100 


60 


L 


5.445 (0.025) 


16^ X 48 


100 


N 


5.5 (0.1) 


24^ X 64 


101 





5.5 (0.05) 


24^ X 64 


100 


M 


5.5 (0.025) 


20^ X 64 


199 


P 


5.5 (0.0125) 


20^ X 64 


199 


G 


5.6 (0.01) 


16^ X 32 


200 


R 


5.6 (0.01) 


24^ X 64 


200 


F 


5.7 (0.01) 


16^ X 32 


49 
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